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The nature of fluorescence intermittency for semiconductor quantum dots (QD) and single 
molecules (SM) is discussed from the viewpoint of Anderson localization. The power law distri- 
bution for the on time is explained as due to the interaction between QD/SM with a random 
environment. In particular, we find that the on-time probability distribution behaves differently 
in the localized and delocalized regimes. They, when properly scaled, are universal for different 
QD/SM systems. The on-time probability distribution function in the delocalized QD/SM regime 
can be approximated by power laws with exponents covering —2 < m < 0. 
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Recent development in nano-fabrication and measure- 
ment has made it possible to probe directly the dynamics 
of a single quantum system and its couphng with local en- 
vironment. For optical measurements, in particular, flu- 
orescence intermittency (FI) occurs in a wide class of sin- 
gle quantum systems including semiconductor nanocrys- 
tals 0, S i, i, and single molecules H S 0. This 
optical phenomenon is experimentally found to be a ro- 
bust and fundamental property of single photon emitters. 
The random switching between an emitting (on) and a 
non-emitting {off) state is characterized by 077.-time and 
o^-timc distribution functions. Surprisingly, the distri- 
bution functions can often be fitted by power laws instead 
of exponentials. The exponents of the power laws vary 
from system to system and they are almost independent 
of temperature. To explain this unusual behavior, Mar- 
cus et. al. 0,13 recently proposed a model based on dif- 
fusion in energy space. An exponent of —1.5 is readily ob- 
tained for one-dimensional diffusion processes. Deviation 
from —1.5 is explained by anomalous diffusion, which has 
its origin from the Cole-Davison dielectric medium /? ^ 1 . 
However, this scenario would predict exponents that are 
temperature dependent. Recent experiments also found 
that exponents ~ —2.0 for the on time exist in many dif- 
ferent single emitters 00,IE0I1|> which are imbedded in 
various media. The temperature-independence of power 
laws with mon/off —1-5 and robustness of rrion = —2.0 
need to be explained. 

In this work, we provide a microscopic Hamiltonian 
model which does not build on the assumption of anoma- 
lous diffusion in energy space and the quadratic form of 
diffusion potential. The purpose of this work is to lay the 
foundation for a quantum-mechanical theory of FI. In a 
classic paper .11.] , Anderson considered a transport pro- 
cess which involves solely quantum-mechanical motion in 
a random lattice. His quantum-mechanical treatment of 
the process leads to the concept of Anderson localiza- 
tion, which has been applied extensively to the problem 
of transport in a random environment. Here, we view 
the blinking as a charge hopping process between the 
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QD/SM and the traps in the environment. We find that 
the FI can be explained in the framework of the Ander- 
son localization. The basic idea is the following: due to 
the coupling to a random environment, the on-site en- 
ergy of a QD/SM is renormalized and becomes a random 
variable. In addition, random configurations of the local 
environment lead to a distribution of decaying (hopping) 
rates for the QD/SM. This is the physical origin of the 
on-time distribution. 

Based on experimental results, we study the following 
simple tight-binding Hamiltonian: 

H = £dd)d + ^ ejc]cj + ^ {d}cj + c]dj (1) 

j 3 

The ground state is defined such that the QD/SM is neu- 
tral. For a QD this corresponds to a filled valence band. 
When an electron is excited by photon to an excited state, 
it relaxes to the lowest excited state on a time scale from 
hundreds of femtoseconds to picoseconds. This lowest 
excited state is responsible for photon emission and hop- 
ping to traps in the matrix. This state could be IS" state 
for QDs and it may involve some surface states as well 
ll^ . The excited electron completes a radiative cy- 
cle on a time scale of the order inverse Rabi frequency: 
^Rabi- This time scale is much smaller than the typi- 
cal bin size, >0.2 ms, used in the experiments. Since 
the excited electronic state is near resonance with the 
ground state plus a photon, we define the on-state as the 
following: it is a dressed electronic state 0| which has 
the same bare energy as the excited electronic state, Ed, 
but with an infinite lifetime when it is decoupled from 
the trap states. Physically, this means that when there 
is no coupling to the environment, the excited electron 
localizes inside the QD/SM and the QD/SM, being in 
the on-state, keeps emitting photons. Such an on-state 
is created by operator . The bare on-site energy £d is 
renormalized due to interaction with trap sites j. The 
coupling constant is Vj. The trap site j is created by 
operator and has a random energy Sj. Here we focus 
on the class of random environments that are formed by 
topologically random network of chemical bonds. Ran- 
domness could arise from a spatially fiuctuating poten- 
tial due to charged impurities and coupling with phonons 
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that arise from deformation of a random lattice. The de- 
tails of randomness are irrelevant to the dynamical re- 
sponse of the system as is known from scaling theory of 
localization. For simplicity, we assume all traps follow 
the same energy distribution function p{s). 

Mathematically, the above definitions are equivalent 
to defining the zeroth order Green's function of the on- 
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state as: G^Vw) — , . , 

being an infinitesimally small and positive constant. Sim- 
ilarly, the zeroth order Green's functions of the traps 



The decay rate of the on- 
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state can be calculated from the full Green's function 
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One obtains the renormalizcd energy and decay rate from 
the pole equation. The results to 0{v^) are: 
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Since each Sj is a random variable, both the renormal- 
izcd energy Ed and decay rate T are random variables. 
Physically, this corresponds to the situation in which en- 
vironment stays at a certain configuration for a period 
of time during which a physical decay rate can be de- 
fined. The ensemble of configurations that the environ- 
ment takes on gives rise to a distribution of decay rates. 
This picture is valid when static disorder dominates or 
when the environment is ergodic. 

To find the decay rate distribution function, /(F), from 
p{sj), one must perform a summation over random vari- 
ables Sj. In addition, due to lattice deformation, a ran- 
dom spatial distribution of traps must also be assumed. 
For simplicity, we assume that all the traps follow the 
same distribution function p(r). We will consider two 
physical regimes, namely the delocalized and the local- 
ized regimes, for the excited electron on the QD/SM ac- 
cording to the framework of Anderson ITlT| . 

The delocalized regime corresponds to most of the ex- 
perimentally studied situations. Essentially, electrons ex- 
cited inside the QD/SM can make real transitions to 



the trap sites in the matrix through the mechanism of 
quantum tunneling, and then become localized. When 
this happens, the QD/SM becomes ionized. The re- 
sults in this regime, which will be shown later, support 
the random resonance picture: on average, the energy 
of a QD/SM, is far off-resonance with trap energies, 
i.e. Ed ^ Sj. The electrons of QD/SM remain local- 
ized {on state) until a random fluctuation of Sj makes it 
on-resonance with the QD/SM. This random resonance 
could be realized through interaction with phonons |17| . 
Since there is finite probability that Sd — Sj ^7/^0"*", 
one cannot perform the summations in the same way 
as in [Tlj . Instead, nontrivial summations of the whole 
expression (|3b|l are required. The distribution function 
/(r) can be calculated by method of Fourier transforma- 
tion: 
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The Fourier transformed probability distribution func- 
tion ip{t) is given by 
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where N is the number of traps and 
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The integral T is obtained under the assumptions made 
above, i.e. identical energy and spatial distributions of 
traps. This integral needs to be done carefully so that 
we can express it in powers of rj. For a real physical sys- 
tem coupled to a disordered environment, rj is small but 
nonetheless not strictly zero. Physically, this is due to 
competing non-radiative relaxation channels which con- 
tribute to a small but finite width of resonance. The value 
of 77 can be calculated in a more sophisticated model, 
however, in our model we take it as a fitting parameter. 
Changing variable to the dimensionless u = — P^tttiS 

and noticing ~ („_„J)2^^2 can be approximated by an 
impulse function centered at u = as 77 is sufficiently 
small, where Ud — 2-Kv{r)'^ \ t\ ^<i^ '^^^ separate the inte- 
gral into three parts: (— 00, Urf — r;/2], {ud~r]/2,Ud + T]/2) 
and [ud + 7y/2, -|-oo) to obtain the expansion: 
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Since we are interested in case where there are a large 
number of traps N ^ 1, hy substituting the expan- 
sion to Eq.© we rewrite the unknown parameter Nij 
as = aQH fill [l8l| . where n is the density of traps 
and ao is a fitting parameter. Absorbing the numerical 
constant in the square bracket of I into oq and neglect- 
ing an unimportant imaginary part as is justified in this 
delocalized regime, we get 



ip{t) = exp 
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The probabihty distribution function /(F) in the delo- 
calized QD/SM regime: 



(7) 



is Lorentzian, with the characteristic decay rate Fq given 
by: 
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The on-time distribution function can be calculated as 
-Pon(i) = /o°° ^rexp(— rt)/(r) which turns out to be: 

= [cos(roi)ci(roi) + sin(Fot)si(Foi)] (9) 

IT 

where ci and si are cosine and sine integrals fl^ respec- 
tively. 

The expression for Pon (t) in has several physical im- 
plications. First of all, it depends only on one parameter 
Fo (this is also true for Pon(i) in the localized regime (|12|l . 
which only depends on the parameter 70). The charac- 
teristic decay rate Fq, is purely due to quantum tunnel- 
ing and depends on the degree of disorder p{ed){v{r)'^) ■ 
This is a generic feature of Anderson localization the- 
ory. Only the energy fluctuation of traps evaluated at 
the bare QD/SM on-site energy, p{ed), comes into the 
final result. This supports the random resonance picture 
discussed above. The one-parameter theory also implies 
that experimental results for different QD/SM systems, 
when plotted in unit of characteristic on-time lifetime 
to = I/Fq, follow a universal distribution Fig^a). In 
the long time limit, t ^ to, Pon{t) ~ t~^. The exponent 
of —2 comes from the peak of the Lorentzian distribution 
for decay rates, which physically indicates the resonant 
scattering condition. This —2 exponent is also robust and 
it does not depend on the properties of matrix or emit- 
ters @. In experiments, the observation time is fixed so 
that one probes only a certain range (3 to 4 orders of 
magnitude in physical time) of the whole probability dis- 
tribution. The location of the window in the universal 
distribution depends on Fq, which is not universal. If 
one fits the experimental data with power laws, then the 
exponents of the power laws will implicitly depend on Fo 
except for the mon = —2 case. From we see that Fo 
depends on properties of the QD/SM and its embedding 



matrix. This is one major difference between our the- 
ory and the diffusion-controlled electron transfer theory 
for the on-time statistics. In the latter, the expo- 
nents of the fitting power laws come from properties of 
the dielectric medium alone, i.e. the f3 parameter of the 
Cole-Davison equation, which depends on temperature. 
As Fo rises with increasing coupling strength between 
the QD/SM and its local environment, power law expo- 
nents observed for SMs would in principle larger than (in 
absolute value) that of well-coated QDs. This is indeed 
consistent with current experimental results 

The quantum mechanical theory of FI allows a weak 
temperature-dependence of the on-time probability den- 
sity. This comes from the fact that the energy fluctuation 
of trap sites, besides the mechanisms discussed above, 
has a temperature component. This dependence is usu- 
ally very weak since the typical energy difference between 
Ed and Ej is on the order of 1 ~ 2 eV which is much 
larger than the magnitude of thermal broadening around 
Sj. Increasing temperature slightly increases p{sd) since 
£d 3* £j and as a result Fo increases. So higher tempera- 
ture would slightly decrease mon towards —2 for the long 
time. 

In the case that the excited QD/SM electron is local- 
ized, then we predict very long on times. This corre- 
sponds to the off-resonance situation, Ed — Sj ^ ?/,Vj. 

One can calculate the probability distribution of F = 

2 

IT (ea-EjV ' using the results of [ll|: 
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W is the width of energy fluctuation for the traps, n is 
the density of traps, while ao ^ 0(1) is a constant with 
a dimension of volume. Note that in the limit 77 — > 0, 
then 7o — > and the probability distribution function of 
decay rates is only nonzero when F = 0, i. e. the QD/SM 
electron is localized. Physically, the smallness of 77 has 
to be compared with the energy difference €d — ej pL5l| . 
As discussed in the delocalized regime, the coupling of 
trap electrons with phonons and non-radiative processes 
etc. can lead to a sufficiently small, but non-zero value 
of rj. The above results are applicable to many topologi- 
cally disordered systems Thus, for a QD/SM in this 
regime, the single emitter experiences very long bright 
periods, which are characterized by 1/70- The on-time 
probability distribution function is obtained from IjlUf) : 
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So in the limit t <^ l/70j the probability distribution 
of the on-time follows a power law with an exponent 
of —1/2. On a much longer time scale, Pon{t) deviates 
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from the power law and shows a stretched-exponential 
tail. This stretched-exponential tail is of purely quantum 
nature and depends on the degree of disorder W/{v{r)). 

Now we discuss the mechanism for the ojf-tmie prob- 
ability distribution. Suppose a charge is on the trap site 
j. The self-energy of Gj{Lo) is Gd{i^)- The decay rate is 
vanishingly small, so the charge is localized. Then on- 
state can only be recovered via quantum tunnelling of 
trapped electrons. Posit) is therefore given by 

-Poff'(t) = / drg{r)^ofiir)exp[-jofi{r)t] (13) 
Jo 

where g(r) is the probability density of an excited elec- 
tron having been trapped at a distance of r from the 
QD/SM and it is proportional to the square of the excited 
state wave function. 7off(f') is the decay rate of a trapped 
electron to the QD/SM. A model based on this picture 
was first proposed by Verberk et.al. Q and experimen- 
tally examined in H. Since the typical size of a QD/SM 
is of the order 30 A in radius, so to a good approxima- 
tion, the QD/SM can be viewed as a shallow impurity 
and one can use effective mass theory. The wave function 
of an excited QD/SM electron behaves as exp(— r/a) at 
large distances, where a = fi^ e^n / m*^e'^ = cik is the ef- 
fective Bohr radius and is related to the ionization en- 
ergy, E'ion, of the excited electron by a = h/ yJ^m^jEi^. 
K is the relative static dielectric constants of the me- 
dia. TO^ is the effective mass of the QD/SM electron. 
Similarly, 7off('") is given by 7off exp(— 2r/6), where 
b = h/ y^2m^Et is the spatial extent of a trapped elec- 
tron's wave function with and Et being the effective 
mass and ionization energy at the trap site. We obtain 

Pos{t) « Ct-™°« with C = ar(l + b/a)/-i'J^b and 

h m*,Eion bl , , 

a y m^Et a k 

We see that the difference in motf for different systems 
comes mainly from the effective mass and the abil- 
ity of a matrix to stabilize the charged QD/SM and the 
ejected electron. In contrast to the diffusion model 0, 
the exponents of Poff(i) depend on the relative static di- 
electric constant k instead of the Cole-Davison /? param- 
eter. This is consistent with recent experimental results 
0, 1^ . In addition, our model allows TOoff to take a differ- 
ent value from mon Q, ^ IE> Q > while the diffusion model 
predicts that they must be the same. 

To conclude, we have proposed a mechanism for FI 
as being a manifestation of Anderson localization. The 
on-time manifold is shown to be generated by different 
realizations of electron delocalization from the QD/SM 
through the mechanism of random resonance. The quan- 
tum theory predicts a universal probability distribution 
function for the on time and shows mon — ~2 is indeed 



a robust result, which corresponds to the long time limit 
of Pon{t). The off state corresponds to a localized elec- 
tron in the environment. The recovery of the on state is 
realized via quantum tunneling. 
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FIG. 1: Probability distribution functions for the on time 
in the delocahzed regime, (a) Universal on-time probability 
distribution function obtained by scaling in a characteristic 
on-state lifetime to = I/Fq. (b) Pon(i) for a QD/SM with a 
lifetime 2ms and a fixed observation window. The exponents 
of the fitting power laws lie between —1.6 and —2. (c) Pon{t) 
of a QD/SM with a shorter lifetime 0.5ms. The exponents of 
the fitting power laws lie between —1.8 and —2. 
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